Accepted for publication in ApJ 

Preprint typeset using I^'T^]X style cmulatcapj v. 04/20/08 



o 
o 

(N 

a 

^3 



6 



> 

O 

o 
o 



X 



TRIGGERED STAR FORMATION IN A DOUBLE SHELL NEAR W51A 
Miju Kang^'^'^, John H. Bieging"^, Craig A. Kulesa''', and Youngung Lee^ 

Accepted for publication in ApJ 

ABSTRACT 

We present Heinrich Hertz Telescope CO observations of the shell structure near the active star- 
forming complex W51A to investigate the process of star formation triggered by the expansion of an 
H II region. The CO observations confirm that dense molecular material has been collected along the 
shell detected in Spzteer IRAC images. The CO distribution shows that the shell is blown out toward a 
lower density region to the northwest. Total hydrogen column density around the shell is high enough 
to form new stars. We find two CO condensations with the same central velocity of 59 km s~^ to the 
east and north along the edge of the IRAC shell. We identify two YSOs in early evolutionary stages 
(Stage O/I) within the densest molecular condensation. From the CO kinematics, we find that the 
H II region is currently expanding with a velocity of 3.4 km s~^, implying that the shell's expansion 
age is ^ 1 Myr. This timescale is in go od agreement with numerical simulations of the expansion of 
the H II region (jHosokawa fc Inutsukal ). We conclude that the star formation on the border of the 
shell is triggered by the expansion of the H II region. 

Subject headings: HII region — infrared: ISM — ISM: bubbles ~ ISM: individual (W5IA) — ISM: 
kinematics and dynamics — stars: formation 



1. INTRODUCTION 

Feedback from massive stars has significant impact 
on the surrounding interstellar medium (ISM) , and may 
even trigger star form ation in the vicinity of their ex- 
panding H II regions (Elme green fc La"dalll977f) . An ex- 
panding H II region sweeps up an ambient interstellar 
medium and accumulates the material between the ion- 
ization front and the shock front. The material collected 
on the border of the H II region will form a shell struc- 
ture and would appear as a ring on the plane of the 
sky. Second-generation stars may form when the com- 
pressed layer becomes gravitationally unstable or the ex- 
panding shell provides enough external pressure to initi- 
ate collapse of a pre-existing molecular clump. There 
are relatively few good ex amples o f this proces s pub- 
lished to date (e.g.JP charvc ng et al.|[2 003: Zavagn o et al.l 
[200l I Watson et al.ll2008i : iPomares et al.ll2009ll . There- 
fore, identifying other cases of shell structures around 
OB stars, together with detailed study of the associated 
interstellar matter, will provide valuable understanding 
of the triggered star formation process. 

High spatial resolution of the Spitzer Space Telescope 
makes it possible to routinely detect these shell struc- 
tures at mid-infrared wavelengths. iChurchwell et al] 
(|2006| ) found 322 shell structures from the Galac- 
tic Legacy Infrared Mid-Plane S urvey Extraordinaire 
(GLIMPSE; IBeniamin et all 120031 ). They reported that 
many of the infrared shells coincide with known H II re- 
gions produced by O and early-B type stars. About 7% of 
the shells in their catalog reveal a multiple morphology. 
The shell structure we present in this paper is a multiple 
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bubble consisting of N102 and N103 in their catalog. 

In this paper, we use infrared and CO observations to 
study the double-shell structures and their interactions 
with the surrounding ISM. Infrared observations from 
Spitzer show the morphology of the shell and reveal the 
YSOs possibly formed by triggering in the surrounding 
ISM. Direct observation of molecular gas around the shell 
is necessary (I) to confirm that ambient molecular mate- 
rial is indeed associated with the expanding H II region, 
(2) to find any clumpy material that could condense to 
form new stars, and (3) to learn the physical properties 
(e.g., density of accumulated molecules) of the ambient 
ISM into which the shell is expanding. These observa- 
tions will enable us to determine whether or not it is 
possible to form stars by the triggering effect of the H II 
region. 

This paper is organized as follows. We introduce 
our CO observations and the data sets used in this 
paper in §[21 In ji l3.Il we describe the morphology of 
the N1G2 and N1G3 region. We identify the ionizing 
star in §[3T2] and any young stellar objects (YSOs) 
around the shell in ii l3.3l In §2] we present the distri- 
bution and kinematics f? l4.ip and physical conditions 
(§53 of the molecular cloud associated with the shell. 
We discuss the triggered star formation on the border 
of the H II region in §[5] We summarize our results in §[6l 

2. OBSERVATIONS AND DATA ANALYSIS 

We have carried out ^^CO and "CO (J=2-I) hue ob- 
servations of the W51 H II region complex with the 10 
meter Heinrich Hertz Telescope (HHT) on Mt. Gra- 
ham, Arizona. The whole map covers a 1.25° x 1.00° 
region centered at {l,h) = (49.375°, -0.2°). In this pa- 
per, we extract a 0.152° x 0.150° region centered at 
(/, h) =(49.696°, -0.148°) which includes the shell struc- 
ture. The i^co and ^^CO (J=2-l) line observations 
were mapped with on-the-fly (OTF) scanning in longi- 
tude at 10" per second with row spacing of 10" in lat- 
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Fig. 1. — A color image of the W51 complex composed of MIPS 
24 iim (red) and IRAC 8.0 fira (green). The contour map shows 21 
cm radio continuum emission. It was made by combining VL A data 
and Effolsborg 100-m telescope data (K oo fc MoonI 1199711 . The 
contour levels are -0.015, 0.015, 0.05, 0.10, 0.3, 0.5, 1.0, 1.5, and 
2.4 Jy beam~^. The white box shows the outline of the shell area 
(Figure [g. 

itude. We used the 1.3mm ALMA band 6 dual polar- 
ization sideband-separating receiver with a 4-6 GHz IF 
band. The receiver was tuned with the ^^CO (J=2-l) line 
at 230.538 GHz in the upper sideband and ^^CO (J=2-l) 
line at 220.399 GHz in the lower sideband. The spec- 
trometers, one for each of the two polarizations and the 
two sidebands, were filter banks with 256 channels of 1 
MHz width and separation. Beam efficiency was 0.85, 
which we adopt for all the data. 

Data for each CO isotopomer were processed with the 
CLASS reduction package (from the University of Greno- 
ble Astrophysics Group), by removing a linear base- 
line and convolving the data to a square grid with 10" 
grid spacing. The intensity scales for the two polariza- 
tions were determined from observations of W51D made 
just before the OTF maps. System temperatures were 
calibrate d by the standard am bient temperature load 
method (jKutner fc Ulichlll981h after every other row of 
the map grid. Further analysis was done with the Miriad 
software package (Sault et al. 1995). Although the data 
observed early in the program with a different receiver 
were noisier than with the ALMA receiver, data qual- 
ity was improved after combining all of the data using 
(l/Tsys)'^ weighting. To further reduce the noise of the 
images without significantly sacrificing angular resolu- 
tion, we convolved the maps with a 16" (FWHM) Gaus- 
sian. This increases the effective resolution from 32" to 
36", but smooths over the 10" OTF raster pattern, im- 
proving the final image quality. The two polarizations 
were averaged for each sideband, yielding images with 
rms noise per pixel and per velocity channel of 0.16 and 
0.07 K-T^ for the ^^CO and ^^CO transition respectively 
The average rms noise level of the ^^CO map is larger 
than ^^CO because several sub-fields of ^^CO map were 
observed just one time with the noisier receiver. 

^^CO (J=l-0) data covering the same region as our 
mapping area were extracted from the Galactic Ring Sur- 
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Fig. 2. — (a) Composite image of IRAC 8.0 fim (red), IRAC 
4.5 fim (green), and IRAC 3.6 fim (blue). Contours show Wester- 
bork Synthes is Radio Telescope(WSRT) 327 MHz radio continuum 
UTavloir et al. 1996). Contour levels are 0.07, 0.12, 0.16, and 0.21 Jy 
beam~^. Scale bar on the left bottom represents a size at the dis- 
tance of 5.7 kpc. (b) Composite image of MIPS 24 /^m (red), IRAC 
8.0 fim (green), and IRAC 4.5 ^m (blue). Outer solid circle repre- 
sents a shell size of 2' radius centered on (Z, b) = (49.698°, -0.162°) 
of N102. Inner dash-dotted circle represen ts N103 with 0.67' radi us 
centered on {l,b) = (49.703°, -0.168°) llChurchwell et al.l 120061 '). 
The white box is the outline enlarged in Figure |3] 

vey (GRS) ^ with spectral resolution of 0.2 km s""'^, an- 
gular resolution of 46" and sampling of 22". The rms 
noise level of the GRS data when smoothed to 1.3 km 
velocity resolution is 0.08 K-T^. The intensities on 
an antenna temperature scale were divided by the main 
beam efficiency of 0.48 to convert to main-beam bright- 
ness temperature. 

The Galactic Legacy Infrared Midplane Survey Ex- 
tradodinaire (GLIMPSE) I survey observed the Galac- 
tic plane (65° < \l\ < 10° for |6| < 1°) with the four 
mid-IR bands (3.6, 4.5, 5.8, and 8.0 fiin) of the In- 
frared Array Camera (IRAC; iFazio et all I2004D on the 
Spitzer Space Telescope. The Midinfrared Imaging Pho- 
tometer for Spitz er Galactic plane Survey (MIPSCAL; 
iCarev et al.l 120051 ) is a legacy program covering the in- 
ner Galactic plane, 65° < |;| < 10° for \b\ < 1°, at 24 
and 70 fj,m wit h the Multiband Im aging Photometer for 
Spitzer (MIPS: iRieke et al.|[200l on the Spitzer Space 
Telescope. We analyzed mosaicked images of four IRAC 
bands and MIPS 24 ^m. The resolutions of IRAC and 
MIPS 24 fim are 1"2 and 2'.'4, respectively. 

We use the GLIMPSE I Catalog which consists of point 
sources that are detected at least twice in one band with 
a S/N > 5. The GLIMPSE I Catalog also includes 
JHKs flux densities from the 2MASS point source cat- 
alog (Skrutskic et al. 2006). MIPS 24 ^m fluxes were 
obtained from running the GLIMPSE pipeline on the 
MIPSGAL BCD mosaics (Meade, private communica- 
tion). We extract 24 ^m point-sources with S/N > 7, 
then bandmerge the 24 fim sources with the GLIMPSE 
Catalog sources using a 2"0 correlation radius. The fl- 
nal source list consists of all 8 bands combined: 2MASS 
JHK,, IRAC (3.6, 4.5, 5.8, and 8.0/im), and MIPS 24 
fj,m. 

3. PROPERTIES OF THE REGION 

3.1. Spitzer Images and Radio Continuum 

Figure [1] shows a two-color composite image of Spitzer 
IRAC 8 fim (green) and MIPS 24 fim (red) data from the 

* http:/ /www. bu.edu/galacticring/ 
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Fig. 3. — Images of the central region at the indicated wave- 
lengths. Candidate ionizing stars are labeled from si to s9 on the 
3.6 iim image. The large circles are centered on N103 and have 
a radius of 0.67'. Two crosses in 4.5 /im image are the center of 
N102 and N103. 

GLIMPSE and MIPSGAL Legacy data bases, for a l°x 
1.25 °area including the W51 star-forming region. The 
contours outhne the radio continuu m emission from th e 
ionized gas and supernova remnant ( Koo &: MoonllTOOTl ) . 
The double shell structure studied in this paper is located 
in the white boxed area of Figure [T] 

Figure [2] shows an enlarged view of this double shell 
structure from the composite image of IRAC and MIPS 
bands. The i nside and outside shells a re N103 and N102, 
respectively (jChurchwell et al.' '2006"). A bright 8 fiin 
shell (solid line in Figure [2l)) encloses bright 24 /im emis- 
sion from the inner shell region. IRAC 8 /im mainly 
shows the emission of polycyclic aromatic hydrocarbons 
(PAHs), which are excited in the photodissociation re- 
gion (PDR). Therefore PAHs emission trace the ioniza- 
tion fronts well. MIPS 24 fim present the emission of the 
small dust grains, which is inside of the shell. We note 
that these shell structures are not visible in any of the 
2MASS near-IR images, despite their prominence in the 
GLIMPSE IRAC data. Average radii of N102 and N103 
are 2.'0 and 0'.67, corresponding to 3.3 pc and 1.1 pc, 
respectively. Here we adopt a kinematic distance to the 
shells as 5.7 kpc, which were deduced from the velocity 
of the associated molecular material, Vlsr = 59 km 
from ou r HHT observa t ions, and the Galactic rotation 
curve of iBrand fc Bliti (fl993l ). 

This shell structure lies in the direction of the ma- 
jor star-forming regions in the Galaxy (Avedisova 2002). 
This apparent association on the plane of the sky does 
not prove that the shell results from the expansion of 
an H II region caused by star formation activity, rather 
than some other process (e.g., supernova remnant, plan- 
etary nebula, etc.). To determine whether our dou- 
ble shell represent is H II region or not, we compile 
the radio continuum observations at various frequen- 
cies: 5'4.s7.qnHz = 0.55 Jy and S^Ann Tjv. = 0.21 Jy 
(IWink et all Il98l . aci ^, = 0.62 Jv (ICondon et all 



[1998D, ^4.8.^ qHz = 0.33 Jy (IGreeorv et al.l ll996fl. ^327 mhz 
= 0.80 Jv (iTayloret al.l 119961). g^sv.nM, = 0.55 Jy 
(jAltenhoff et al.lll979n . 













■ = 10 ^ 










s2 sg^"^/^ 



















































a) - 



0.0 



0.5 



1.0 

H-K, 



1.5 



3.0 



2.0 



1.5 



1.0 



0.5 



0.0 



-0.5 



Av = 20 
t 



A A 



s8 „ 



b) 



-0.5 



0.0 0.5 

[5.8]-[8.0] 



1.0 



Fig. 4. — (a) JHKs color-color plot. The solid and dashed 
curves represent the locus of main-sequence and giant stars 
IIBessell fc Bret j 119881). The parallel solid and dashed lines, for 
06-8V (Koornnccf 1983) and MOV jBesscll & Brett 1988), are rcd- 
dcni ng vectors with color e xcess ratios of 2.07 (Kim et al. 2007) and 
1.9 llOkumura et al.l 120001) . respectively. Reddening vector shows 
Av=10 mag. (b) IRAC [3.6] - [4.5] vs. [5.8] - [8.0] color-color 
plot. Gray dots are all sources detected at each band through the 
line of sight in the shell region. Open diamonds are stars within 
the radius of N103 seen in Figure |3] Filled diamond represents a 
candidate ionizing star. YSO candidates are marked as circles for 
Stage O/I, triangles for Stage II, and squares for ambiguous sources. 
Re ddening vector for Ay =20 mag based on the extinction laws of 
llndebetouw et al.i (i2005l ) is shown as a filled arrow. 



iTavlor et al.l (|1996l ) concluded that most sources de- 
tected by the Westerbork Synthesis Radio Telescope 
(WSRT) at 327 MHz and in the IRAS Point Source Cat- 
alog are thermal radio sources associated with compact 
H II regions. Their W1921-fl442 is coincident with the 
shell structure as seen in Figure [2] ( a) . IRAS fluxes of 
the source are 1.0, 24.9, 225.8, and 29 2.5 Jy at 12, 25, 
60, and 100 /xm. ITavlor et all (|1996D selected a ratio 
of S'eo^^m to 5*327 MHz = 200 as the boundary separat- 
ing thermal and non-thermal radio sources. The ratio of 
5*60^111 to S'327MHz of W1921+1442 is greater than 200, 
consistent with thermal emission. The IRAS flux ra- 
tio also hes in the region 0.5 < log(S'6o/5'25) < 1.7 and 
0.5 < log(525/<S'i2) < 1.4 occupied by compact H II re- 
gions. 

We use equation (1) of 'Si mpson fc Rubinl (jl990l ) to es- 
timate from the radio-continuum flux the number of ion- 
izing Lyman continuum photons A^Lyc emitted by the 
exciting star. We assume a fraction of helium recom- 
bination photons to excited states ji = 0.65, an elec- 
tron temperature = 10000 K and a distance D = 5.7 
kpc. The integrated flux density at 4.85 GHz is 0.33 
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TABLE 1 

Model parameters for YSO candidates^ 
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The 0.152° X 0.150°rcgion centered on (l,b) =(49.696°, -0.148°) 
^ The value for the best-fitting SED. A good fit(x^ /^data ^ 4) can have a value as high as 16 or more. 

^ The uncertainties on the luminosities, masses, mass accretion rates, and extinctions are calculated as the weighted standard 
deviation of the luminosities, masses, mass accretion rates, and extinctions of all the acceptable YSO models. 



Jy (jGregorv et al.|[T99()f) . indicating an ionizing flux of 
log(A^Lyc) ~ 47.96 photons per second. This corre- 
sponds to a star of spectral type 09V ("Martin s et al 
2005). T he value 5'4.87 5GHz = 0.55 Jy (AltcnhofT^eLal 
1979tlWink et al.| [T98a gives log(iVLyc) 48.18 photons 



per second, which implies on 08.5V spectral type. The 
ionizing flux may be a lower limit if there is significant 
dust absorption in the ionized gas. 

3.2. Identifying Ionizing Star 

Given that the N102 and N103 structures are shells 
swept up by expanding H II regions, we try to identify 
the ionizing stars using NIR photometry. N102 (outer 
shell) shows a nearly complete circular shape except for 
the missing region to the north-west which suggests an 
eruption of the bubble. Ionizing star(s) might be located 
near the center of the shell where the inner shell (N103) 
is located. Therefore, we consider stars detected at all 
2MASS bands within the inner shell as possible candi- 
dates for ionizing star(s). Figure [3] shows the candidate 
ionizing stars, labeled from si to s9, at different IRAC 
and MIPS bands. 

The most likely candidate to be the ionizing source 
was selected as follows: Using photometry data from the 
2MASS and Spitzer catalogs, we placed sources in Fig- 
ure [3] on color-color diagrams (see Figure [4|) . Extinction 
is clearly severe toward the Sagittarius spiral arm tan- 
gency at 45° < I < 51° including W51, since t here is 
no detectable increase in the star count (iB cniami n et al.l 
[lOOl iChurchwell et al.' [200l . lOkumura et all (|2000f ) 
and iKim et al. ( 2007 ) estimated a color excess ratio 
Ej_h/Eh-k, of 1.9 and 2.07 toward W51A and W51B, 
respectively. In the J — H vs. H — Kg color-color dia- 
gram, we assumed a color excess ratio of 2.07. Although 
the color excess ratio is measured from the W51B region, 
the value 2.07 traces the sources in the shell region bet- 



ter than 1.9 measured in W51A (see Figure |4ji). This 
excess ratio is at the upper range of 1.73 ib 0.40 found 
in the Galactic reddening law stud y of llndebetouw et al.l 
()2005f ). Following the criterion bv lKim et all ()2007D . we 
rejected sources with J — Ks < 2.2 as foreground sources. 
Notice how this includes source s4, which despite being 
located near the center of the N102 structure, presents 
negligible extinction a nd has a color c o nsiste nt with the 
locus of giant stars bv lBessell fc BrettI ()1988D : moreover, 
the source is easily detectable in optical plates from the 
Digitized Sky Survey (DSS), which would exclude it from 
being a deeply embedded source. Foreground sources si 
and s3 are also marginally detected in the DSS optical 
band. Sources s2, s6 and s7, have J — Kg > 2.2 showing 
significant reddening, but they are located in the region 
of the diagram corresponding to the extension of the cool 
Giant locus towards highly reddened values, and thus we 
also discard them as candidate ionizing sources. The re- 
maining source, s8, is located in a region of the diagram 
likely to be occupied by deeply embedded O type stars. 

Candidate ionizing stars detected in all IRAC bands 
are located in the group of main-sequence and giant stars 
around (0,0) in the Spitzer [3.6] - [4.5] vs. [5.8] - [8.0] 
color-color diagram in Figure H] ( 6) , indicating that they 
are not YSOs. 

3.3. YSO Candidates 

In order to determine whether the expanding H II re- 
gion is able to trigger the star formation, we investigate 
if there are newly formed stars in the vicinity of the shell 
structure. To find YSO candidates associated with the 
shell structure we use the 2MASS and GLIMPSE point 
source catalogs. There are 2189 sources in 0.152° x 0.150° 
region around N102. We identity a t otal of 19 YSOs 
near the shells using the SED fitter of 'Robitaille et al 
(j2007.) following the procedures described in .Povich et al 
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Fig. 5. — Top: Average I^CO and ^^CO (J=2-l) spectra seen 
in Figure [2] Bottom (a) The ^^CO (J=2-l) intensity integrated 
between 40 and 50 km s-^. (&) The ^^CO (J=2-l) intensity inte- 
grated between 56 and 67 km s~^ . The gray scale is from to 44 K 
km s~^. Contour levels are 5, 10, 15, 20, 25, 30, and 35 K km s~^. 
Circles show the shell size of N102. The dash-dotted rectangles 
labeled 1 and 2 show the fields of clumps 1 and 2 in Figure |6] 



([2009f l. The SED fitter includes a grid of 200,000 YSO 
model spectral energy distributions (SEDs) and finds a 
model fit using a ^'^-minimization. For SED fitting, we 
use the sources detected in at least four bands {Ndata > 
4) of the 8 IR bands, consisting of the 4 IRAC bands 
and the 3 2MASS bands for sources with firm 2MASS 
identifications in the GLIMPSE Point Source Catalog, 
and the MIPS 24 ^m flux from the MIPSGAL survey if 
detected (see Section[2|). First, to remove stellar sources, 
we fit stellar photosphere SEDs to all the sources with 
Ndata > 4. Sources with x^/^data < 4, indicating a 
good stellar fit, are removed. Because the fittin g tool 
take into accounts extinction (Inde betouw et all 12005) . 
even highly reddened stars should be removed by this 
step. Asymptotic giant branch (AGB) stars are removed 
by fitting AGB star SED templates included in the SED 
fitting tool. 

Next, we fit the sources which are p oorly fit by stel- 
lar photospheres to the YSO models of iRobitaille et aT] 
(|2b06i) . We allow the distance range to be from 5 to 9 kpc 
and the interstellar extinction from to 60 mag in V band 
for fitting YSOs. In the same manner as fitting stellar 
photospheres, we consider sources with /Ndata < 4 to 
be good fits. We discard sources not detected at 24 /j,m 
which show 8 /im flux excesses above an extrapolation 
of the 3 shorter wavelength IRAC bands. These sources 
are usually contaminated at 8 /im by a noise peak or a 
diffuse background feature. 

Each YSO stage is defined by a disk mass Mdisk 
and an envelope accretion rate Menv provided by the 
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Fig. 6.— I^CO (J=2-l) (thick solid line), ^'^CO {J=2-l) (thin 
solid line) and ^^CO {J=l-0) (dash-dotted line) spectra of clumps 
l{left) and 2 (right) seen in the dash-dotted boxes in Figure [5] The 
spectra have 22" spacing. 

SED fitting algorithm: Sources are classified as Stage 
O/I with M^nv/M^ > 10^6 yj.-i^ g^age II YSOs with 



M^nv/M^ < 10-^ yr^i and M^isk/M^ > IQ-^ , and Stage 

III with Menv/Mvc < 10"^ yr"^ and Mdisk/A^* < 10"^ 
(jRobitaille et al.ll2006f ). We determine the evolutionary 
stage of each source using the relative probability distri- 
bution for the Stages of all the good-fit models. Note that 
we use the term "Stage" to distinguish the evolution- 
ary state of the theoretical model from the term "Class" 
which derives from the observed SED (Evans et al. 
The good-fit models of each source are defined by 



X 



A. mm 



< 2iVdata, 



(1) 



where Xmin ^^i^ goodness-of-fit parameter for the best- 
fit model. The relative probability of each good-fit model 
is estimated according to 



-(x^-xL„)/2 



(2) 



and is normalized. After a probability distribution for 
the evolutionary Stage of each source is constructed from 
the Stages of all the good- fit models, the most prob- 
able Stage of each source is determined by requiring 
SP(Stage) > 0.67. If this condition is not satisfied, then 
the Stage of the source is considered as "Ambiguous". 
Model parameters of all YSOs are listed in Table 1. 

FigureH) shows IRAC [3.6] - [4.5] vs. [5.8] - [8.0] color- 
color diagram of sources within the shell region. Note 
that we show only sources that are detected in all IRAC 
bands. The gray dots represent foreground or back- 
ground stars. Candidate ionizing sources seen in Fig- 
ure [3] are marked with black diamonds. Filled diamond 
represents the most probable ionizing star. YSOs with 
IR excess are marked as circles for Stage O/I, triangles 
for Stage II, and squares for Ambiguous. In the IRAC 
color-color diagram, YSO candidates are clearly sepa- 
rated from stellar sources around (0, 0) where candidate 
ionizing stars are located. In the §[Sl we will discuss YSO 
candidates associated with shell. 

4. MOLECULAR CLOUDS 
4.1. Spatial Distribution and Kinematics 
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Fig. 7. — The ^-^CO (J=2-l) channel maps by 1 km s ^ interval. Black contour shows the shell structure of IRAC 8 fim. Contour level 
is 50 MJy sr~^. In last channel map, red contours are MIPS 24 fiui. Contour levels are 50, 100, 200, and 400 MJy sr"!. The color scale 
bar is presented in units of K-Tmb . 



To determine some physical properties of the sheU 
structm-e, we used the ^^CO (J=2-l), ^^CO (J=2-l), and 
^^CO (J=l-0) maps described in Section 2. The first 
step is to determine the velocity range of the molecular 
components associated with the shell structure. Figure 
El (top) shows the velocity profile of each line averaged 
over the whole region shown in Figure [2] There are two 
bright components around 40 - 50 km and 56 - 67 
km s~^ in the ^^CO (J=2-l) line. We present the inte- 
grated intensity maps of ^^CO (J=2-l) in two different 
velocity ranges in bottom of Figure [51 

In the ^^CO (J=2-l) map integrated from 40 to 50 km 
s~^ (Figure llh) , a strong component is located on the 
north-east side of the shell. Some bright emission is seen 
within the radius of N102. The ^^CO (J=2-l) intensity 
map integrated from 56 to 67 km (Figure (5?*) shows 
that some clumpy structures surround the ionized region 
seen in the IRAC shell. We label two of these as clumps 
1 and 2. The ^^CO , "CO (J=2-l) and ^^CO (J=l-0) 
spectra around clump 1 and 2 are presented in Figure [H 
Peak intensity of ^^CO (J=2-l) appears around 59 km 
s~^ in both clumps. The ^^CO (J=2-l) line between 56 



and 62 km s ^ of clump 1 is very bright and distinguished 
from the other components in different velocity ranges. 
The ratio, i^/is^^^-i, of ^^CO (J=2-l) and "CO (J=2- 
1) at the peak velocity of "CO (J=2-l) is 1.85 at the 
position marked as an asterisk in Figure [6] {left). In all 
spectra of clump 1 and 2, the main beam temperature 
of GRS ^^CO (J=l-0) is brighter than the main beam 
temperature of HHT ^'^CO (J=2-l) implying the non- 
LTE state of shell region (§4.2). 

Figure [7] shows "'^^CO (J=2-l) emission associated with 
the shell between 56 and 67 km s~^ by 1 km s~^ steps. 
The molecular and IRAC shells show good correspon- 
dence between 58 and 63 km s~^ in the north-east quad- 
rant where the shell is complete. To the north-west the 
shell is broken in both the IRAC and molecular line 
images. To the south the correspondence is less clear 
because of the presence of additional molecular gas ex- 
tending off the map in Figure [71 MIPS 24 /^m emission 
represented in the last channel, 67 km s^^, fills the in- 
side of the shell. The distribution of the CO molecular 
cloud agrees with the broken shell structure seen in IRAC 
bands. A molecular clump in the south forms the bright 
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Fig. 8. — ^■^CO spatial and position- velocity diagrams around 
the shell region, (a) is the ^^CO (J=2-l) intensity map integrated 
between 55 and 68 km s~^. Contour levels are 5, 10, 15, 20, 25, 30, 
35, 40 and 45 K km s~^. The cr oss symbol represents th e center 
position of the N102 shell listed bv lChurchwell et al.l l(200^ . White 
circle represents the shell size of 2 arcmin radius. The oblique line 
AB shows the section for the PV map presented in (d). (ft) and (c) 
are position- velocity maps in latitude and longitude across the solid 
lines in (a). The contour levels are {1, 2, 4}xl.6 K-Tmb (Scr). (d) 
A position-velocity diagram along AB. The x-axis is the distance 
from the center of the shell. Contour level is 1.6 K-Tmb- The 
velocity variation expected from equation (1), adopting vq = +3.4 



km s ^, R 
dashed line. 



2 arcmin, and Va 



+59 km s , is shown by the 



rim of the IRAC shell. 

Figure [8] shows the spatial and kinematic structure 
around the shell in ^^CO (J=2-l). Figure [51a is the 
^^CO {3—2-1} intensity integrated between 55 and 68 
km s~^. Figures [8j) and [He are position- velocity maps 
in latitude and longitude across the solid lines in Figure 
[SJa, respectively. Figure [SJi is the position- velocity map 
along a line AB in the intensity map (Figure [Sh)- The 
cross represents the center position of the N102 shell. Al- 
though most CO emission is distributed broadly from 55 
to 67 km s^^, the shell structure is visible in the velocity 
range Vlsr — 56 to 62 km s"-^. 

For a simple shell structure, the expected velocity 
along the line of sight can be expressed as 



v{X) 



(3) 



where vq is the expansion velocity, R is the radius of 
the shell (2') estimated from the IRAC images of N102, 
X the projected distance from the center of the shell, 



and 



the systematic LSR velocity of the whole sys- 



tem. We assume Vsys — -1-59 km s since all clumpy 
structures on the border of the shell show the peak at 
Vsys — +59. The expansion velocity of the shell is deter- 
mined to be Vq ~ 3.4 km by least-squares fitting for 
the mean velocity within the shell radius. The expected 



velocity from equation ([3]) is presented as a dashed line. 
Although molecular material around the center of N102 
has disappeared, especially toward north-west, molecular 
gas far from the central star still remains. The dynami- 
cal time-scale of N102 can be estimated by dividing the 
size, 2', by the expansion velocity vq. Assuming a Vsys 
between 58 and 59 km s^^ for N102, we estimate an age 
of 0.71 to 0.96 Myr for the derived expansion velocity, 
which varies from 3.4 km s^^ to 4.5 km s^^. 

4.2. Physical Conditions 

As mentioned in ii|4.1l the observed ratio of ^■^CO (J=2- 
1) and "CO (J=l-0), "i?2-i/i-o~0.5, shows that the 
shell region is not dense enough to satisfy the LTE con- 
dition. Considering the location of the shell region near 
the active star-forming complex, W51A, the kinetic tem- 
perature must be higher than 12K, derived from the 
maximum main beam temperature of ^^CO (J=2-l). 
We therefore need to use a non-LTE statistical equilib- 
rium treatment of the CO molecular excitation to study 
physical properties of the molecular cloud. We use the 
escape probabil ity radiative t ransfe r and photodissoci- 
ation model of iKulesa et al.l (|2005| ). also described in 
iPovich et al.l ()2009l ). to calculate grids of CO level pop- 
ulations for wide ranges of volume densities (10^ — 10'' 
cm^"^) and temperatures (5 — 300 K) assuming detailed 
balance and steady state. From these model grids, the 
total CO column density at each observed pixel is com- 
puted from the peak temperature, line widths and inte- 
grated intensities of the observed CO lines. Assuming 
that the CO heating is dominated by photon processes 
(e.g. the photoelectric heating of dust), a coarse estimate 
of the incident radiation field can be made for each point 
in the map. The photodissociation model is then applied 
to estimate a total hydrogen column density from the 
CO data. This calculation is based on the CO and Ho 



photodissoc iation treatments of Ivan Dishoeck fc Blackl 
ifTgSSi) and iBlack fc van Dishoeckl (|1987f) . respectivelv. 



using a tot al interstellar carbo n abundance of C/H 
2.4 X 10"'' (jCardelh et al.lll996f ). Because the CO abun- 
dance is a strong function of column density, blind appli- 
cation of a uniform CO "dark cloud" abundance (10~^) 
can lead to a gross underestimate of the total hydrogen 
column and gas mass. The visual extinction Ay, is cal- 
culated h2HLi^I^lusingiV(^f) Mi/ = 1.9 X lO^i cm-2 
mag-i ()Bohhn et al.l[T97l . 

To determine how much molecular material is swept 
by the expanding H II region, we calculate the column 
density of CO between 56 and 62 km s~^ centered on 59 
km . This range corresponds to the velocities that 
the clumps 1 and 2 on the shell occupy on the PV map 
(Figure [3). 

We apply this analysis to the combination of 
13C0 (J=l-0) data from GRS, and the ^^CO (J=2-l) 
and i^CO (J=2-l) data from the HHT. We convolve aU 
spectral line maps to the GRS resolution of 46". From 
the escape probability model, we estimate a total mass 
of 3.2 X 10"'^ Mq distributed from 56 and 62 km s~^ in 
the 0.152° X 0.150° region shown in Figure[2] 

To avoid degrading the resolution of our HHT maps, 
we also apply the non-LTE model to the ^^CO (J=2-l) 
and ^'^CO (J=2-l) data. In this case, the total mass in 
the same velocity range and region is computed to be 
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Fig. 9. — Hydrogen column density contours superimposed on the 
IRAC 8.0 /im image (inverted grayscale). The hydrogen column 
density is integrated over the velocity range between 56 and 62 
km s-^ Contour levels are 40, 50, 60, 70, 80 and 90% of the peak 
hydrogen column density of 2.43 X 10^ cm~^. YSO candidates arc 
marked as circles for Stage O/I, triangles for Stage II, and squares 
for ambiguous sources. The dashed box shows the region in Figure 

[To] 

2.9 X 10"' Mq. The derived total mass from two lines 
of i^CO and I^CO (J=2-l) is - 10% smaller than that 
derived from three lines of ^^CO and ^^CO (J=2-l) and 
^■^CO (J=l-0). Since two results are similar, we present 
the total hydrogen column density N(H) distribution de- 
rived from the HHT ^^CO and ^^CO (J=2-l) with 32" 
spatial resolution in Figure [9] 

5. DISCUSSION 

Using the CO data and the identified YSO candidates, 
we investigate whether our shell region are consistent 
with the predictions from the star formation triggered by 
the expansion of an H II region. Figure [5] shows the re- 
sulting distribution of the total hydrogen column density 
N{H) with the shell structure on the IRAC 8 fim image. 
YSO candidates are indicated by circles for Stage O/I, tri- 
angles for Stage II, and squares for ambiguous sources. 
Total hydrogen column density is derived from 56 to 62 
km s~^, which is the velocity range affected by the ex- 
pansion of the H II region. Most dense condensations are 
distributed along the 8 /zm shell. The distribution of the 
hydrogen in the inner region of the large shell is coinci- 
dent with the distribution of the MIPS 24 fim emission 
encompassed by the IRAC 8 /im shell. The distribu- 
tion of gas condensations along the 8 ^m shell, including 
clump 1 and 2 (see Figure [5]), supports that the molecu- 
lar material collects on the boundary of the shell during 
the expansion of the H II region. 

Accordi ng to the "c ollect and collapse" model 
(|Elmcgrc en fc Ladal Il977t ). a thin layer of compressed 
neutral material forms between the ionization front and 
the shock front as the H II region expands. This layer 
may b ecome gravitationally un stable, fragment and form 
stars. I Whit worth et al.) (|1994f ) have shown that the shell 
collected by the expansion of an H II region will fragment 




49-75 49-74 49.73 
1 (deg) 



Fig. 10. — Ay contours superimposed on the composite image of 
MIPS 24 fim (red), IRAC 8.0 ^lm (green), and IRAC 4.5 /xm (blue) 
of the dashed box region seen in Figure [9] Contour levels are {5, 
6, 7, 8, 9, 10, 11, 12} mag. The candidate YSOs 14 and 16 are 
marked with circle (Stage O/I) and square (Ambiguous). The SED 
model fits to the two sources are plotted on the right. The dashed 
curve shows the reddened photosphere of the central star. 

when its column density reaches 6 x 10^^ cm~^. Based 
on our escape probability model using the ^^CO and 
^■^CO (J=2-l) lines, the hydrogen column density of the 
most shell region is already greater than 1.0 x 10^^ cm~^ 
(the lowest contour in Figure [9] ; 40% of the peak hy- 
drogen column density). Note that this N{H) are a 
lower limit because the CO molecule may be depleted 
by freezeout onto dust gra ins. Therefore by the criterion 
of Whitwo rth et all ([1991 . we conclude that the column 
density of the shell region is sufficient to form a new star 
by the collect and collapse process. 

Dynamical expan sion of the H I I regi on has been sim- 
ulated by Hosokaw a fc Inutsukal ()2006D . They describe 
the results of several models with different central stars 
(12 - 101 Mq) and ambient densities (10^ - 10^ cm'^). 
We estimate that the N102 region has an 08.5V ioniz- 
ing star with the ambient number density of about 10"^ 
cm^"^, which we derive from N{H) of 1.0 x 10^^ cm~^ and 
a line-of-sight depth of 2 pc, corresponding to the pro- 
jected size of Clump 1. We compare the physical proper- 
ties of N102 with the model S19 of lHosokawa fc Inutsukal 
(|2006f ) : central star mass of 19 Mq with an ambient den- 
sity of 10'^ cm^"^. In this model, a star produces a shell of 
radius 3.5 pc after 1 Myr and an unstable region appears 
at about 0.5 Myr. For the N102 shell, the current radius 
of N102 is 3.3 pc at 5.7 kpc distance. Assuming a Vsys of 
59 km s^^ for N102, we estimate an expansion velocity 
of 3.4 km s^^ and an expansion time of 0.96 Myr. The 
shell size and the time scale are in excell ent agreement 
with those of model S19 by Hoso kawa fc In utsuka (2006). 
Therefore, we conclude that triggered star-formation on 
the border of N102 is possible by the expansion of the 
H II region. 

Given that the accumulated molecular material is suf- 
ficient to form stars, we find two YSOs (14 and 16) as- 
sociated with the Clump 1 that is the densest region 
along the shell structure. The hydrogen column density 
at the peak of clump 1 is 2.4 x 10^^ cm~^, correspond- 
ing to the visual extinction of 14 mag. The mass of 
clump 1 shown in Figure [TUl is 2.3 x 10'^ Mq based on the 
escape probability modeling using ^^CO and ^'^CO {3=2- 
1) lines. This mass are lower limits because the CO 
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molecule may be significantly depleted by freezeout onto 
dust grains in the core of the cloud. 

Figure [TOl shows the Ay distribution of clump 1 on the 
MIPS 24 fim (red), IRAC 8.0 fim (green) and 4.5 fim 
(blue) composite images. YSO 14 and 16 are marked as 
a circle and square, respectively. YSO 14 and 16, classi- 
fied as Stage O/I and Ambiguous, lie close to the central 
region of clump 1. MIPS 24 iim emission fluxes of 94 and 
15 mJy are detected from YSO 14 and 16. We also show 
the SED model fits to these source in Figure [TOl YSO 
14 is an intermediate YSO of ~ 6 Mq with an accre- 
tion rate of ~ 8 x 10~^Mq yr~^. YSO 16 has a mass of 
~ 3 M0 with an accretion rate of ~ 1.6 x 10~^Mq yr~^ . 
While YSO 16 is classified as an ambiguous stage, a early 
stage is indicated because the good-fit models include 
62% (Stage O/I), 20% (Stage II) and 18% (stage III). 

If considering the total column density, age and size 
of the H II region predicted from the collect and collapse 
picture, star formation should be possible everywhere the 
material is collected on the border of the shell. As seen 
in Figure [U although the CO emission is collected on 
the periphery along the IRAC shell, we cannot be sure 
whether all the YSOs seen in Figure[5]are associated with 
the material collected by the expansion of the H II region. 
At least, YSO 14 and 16 seem to be associated with the 
collected material directly. We conclude that those two 
YSOs have formed probably by the fragmentation of the 
collected material swept up by the expansion of the H II 
region. 

It may be that the exciting star s8 which we have iden- 
tified was itself formed by such a triggering process. The 
molecular gas in the shell has a very similar velocity to 
that of the gas associated with the W51A H II regions, 
so these clouds may be physically associated. If so, wc 
may speculate that expansion of the ionized regions in 
W51A triggered the collapse leading to the formation of 
0-star s8 in a kind of sequential star formation. 



6. CONCLUSIONS 

We have studied triggered star formation near the mul- 
tiple shells of N102 and N103 (Churchwcll ct al."200l) 
using Spitzer IR and HHT CO observations. We have 
identified the candidate ionizing star (possibly 08.5V) 
of the double shell. The CO observations confirm that 
dense molecular material has been collected along the 
shell which has been detected in the Spitzer IRAC im- 
ages. The CO distribution shows that the shell is blown 
out toward a lower density region to the northwest. We 
find two clumpy CO condensations with the same cen- 
tral velocity of 59 km s~^ to the east and north along the 
edge of the IRAC shell. Total hydrogen column density 
around the shell is high enough to form new stars and 
we identify two YSOs in early stages (Stage O/I) within 
the densest molecular clump 1. Using the CO position- 
velocity map, we find that the H II region is currently 
expanding with a velocity of 3.4 km s~^, suggesting the 
shell's expansion time of ~ 1 Myr. This timescale is in 
good agreement with numerical simulations of the expan- 
sion of the H II region (Hosokawa & Inutsuka 2006) . We 
conclude that the star formation on the border of N102 
is triggered by the expansion of the H II region. 
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